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PACS numbers: 75.10.Jm; 75.50.Ee; 75.40.-s 

We describe the four kinds of behavior found in two- 
dimensional isotropic quantum antiferromagnets. Two of 
them display long range order at T = 0: the Neel state and 
the Valence Bond Crystal. The last two are Spin-Liquids. 
Properties of these different states are shortly described and 
open questions are underlined. 



I. INTRODUCTION 

Exact diagonalizations are restricted to small spin sys- 
tems (up to about TV = 36 spins with the present day 
computers). Some problems as the study of critical 
regimes are out of reach of such an approach, on the con- 
trary we think that this approach might help in answering 
other subtle questions as for example symmetry breaking, 
order out of disorder, qualitative nature of the ground- 
state and of the first excitations of two-dimensional (2d) 
spin- 1/2 magnets. 

We have used it extensively to study the competition 
of long range order (LRO) versus " quantum disorder" at 
T ~ Q. Up to now, by looking to the spectra of model 
spin systems, we have been able to characterize four kinds 
of behavior: Neel LRO (coUinear or not), dimer LRO and 
two different Spin-Liquids. In this paper we will give a 
brief picture of these four generic cases of 2d magnetism, 
underline the more efficient numerical criteria to deter- 
mine the nature of a given unknown system and try to ex- 
tract from the present numerical knowledge what are the 
conditions for the appearance of the four different phases. 
The key point of our approach is the analysis of the sym- 
metries and finite size scaling of the low lying eigen- levels 
of a given hamiltonian on small samples with periodic (or 
eventually twisted) boundary conditions. The evolution 
of the spectrum with the size and shapes of the sample 
has to be thoroughly understood before conclusion: it 
may be a lengthy and sometimes subtle process but it is 
usually quite valuable. Some examples will be quoted in 
the following sections. 



II. NEEL LONG RANGE ORDER 

There is no mathematical proof that Neel long range 
order exists at T = in 2d spin-1/2 systems. But there 
is a very large consensus that it should be the case for 



the Heisenberg nearest neighbor hamiltonian on bipartite 
lattices m. Exact diagonalizations and/or Monte-Carlo 
calculations were performed for the Heisenberg model on 
the square lattice |0|, on the triangular lattice [0,0 and 
on the honeycomb lattice [B . On each of these lattices 
the Heisenberg model displays long range order with a 
sublattice magnetization reduced to 60% of the classical 
value on the square lattice (~ 40% on the triangular lat- 
tice [|,|, 44% on the honeycomb lattice [U|). These 
values of the order parameter might be in error by a 
few percent: in these antiferromagnets with linear Gold- 
stone modes, the scaling law for the magnetization of 
samples of linear size L ^ N^'''- is an expansion in 1/L 
(i.e. 1/A^^/^ in 2d) 0^- For the sizes encountered in ex- 
act diagonalizations, the asymptotic l/N^^"^ law is never 
reached ||,|l^ and the extrapolation to the thermody- 
namic limit remains uncertain: nevertheless the agree- 
ment between the scaling behaviors of diagonalization 
results and spin wave results gives a strong support to 
the above-mentioned conclusions. 

A stronger confirmation of Neel LRO is obtained by 
the scaling of the spectrum itself, which is an illuminat- 
ing illustration of the mechanism of symmetry breaking. 
In order to understand in a simple way how a continuous 
symmetry breaking manifests itself on finite size spectra, 
let us consider a piece of solid. Our experience of the 
macroscopic solid is the following: it breaks the transla- 
tion symmetry and its center of mass (c.o.m.) is localized 
in space. On the other hand, we can sort its degrees of 
freedom in two parts: first the c.o.m. coordinates and 
the internal variables. The hamiltonian of such a prob- 
lem breaks in two commuting parts: the c.o.m. kinetic 
energy and the hamiltonian of the internal degrees of 
freedom ( in a solid, in a first approximation this second 
part is a phonon hamiltonian). For a finite size sample 
of N atoms, with mass m, this reads: 



Htnt — 



2mN 



^phonons- 



(1) 



The total hamiltonian docs not break translation symme- 
try, the total momentum P is a constant of the motion 
and its eigenvalues good quantum numbers. Due to the 
commutation property of the sub-parts of the hamilto- 
nian, the eigenstates are direct products of eigenstates 
of the c.o.m. kinetic energy and phonons eigenstates. 
From these very simple observations it follows that the 
spectrum of Eq. 1 as a function of the eigenvalues of the 



square of the momentum has the appearance described in 
Fig. 1. A straight hne is drawn across the eigen-energies 
of the c.o.m. hamiltonian with phonon (respectively 
1 phonon...). Let us first concentrate on the vacuum of 
phonons: the associated eigen-levels collapse to the abso- 
lute ground-state as 1/A'' when N goes to infinity. Let us 
now consider the eigenlevels immediately above this fam- 
ily: they correspond to the softest phonon excitation and 
are distant from the vacuum family by the energy of this 
phonon, that is an energy 0{\/ L) R When the size goes 
to infinity, if the dimension of the solid is equal or larger 
than two, the levels forming the vacuum of phonons col- 
lapse to the absolute ground-state faster than the soft- 
est phonon! This allows to form a wave packet out of 
these levels and localize the c.o.m. with a cost in energy 
lower than the softest phonon: the solid phase with its 
specific macroscopic properties (symmetry breaking and 
stiffness) can exist in 2d (at T = 0), whereas it cannot 
exist quantum mechanically in Id. 

Let us now turn to the Neel ordered magnet. For 
simplicity we will first consider a coUinear magnet: the 
macroscopic state of the magnet at T = is defined by 
the Euler angles of the order parameter. They are the 
equivalent of the c.o.m. coordinates in the solid prob- 
lem. The conjugate variables of these Euler angles are 
the components of the total spin and their free preces- 
sion is described by the Hamiltonian: 



(2) 



2xN 



This is the equation of motion of a rigid rotor. The mo- 
ment of inertia of this rotor is written here in terms of 
the intensive homogeneous spin susceptibility x, -^ is the 
number of spins. If the system has LRO and rigidity, 
the motion of the internal variables is governed by an 
hamiltonian of (antiferromagnetic ) magnons and the to- 
tal hamiltonian reads: 

§2 

^tot — 7^ Tw" '^ 'T'maqnons ' \'^) 

2xN 

When Htot is SU{2) invariant, S^ is a constant of motion 
and its eigenvalues S{S + 1) are good quantum num- 
bers. It is thus interesting to display the spectrum of 
this ordered antiferromagnet as a function of S'(5' -I- 1). 
Its generic appearance is shown in Fig. 2. One recog- 
nizes the (Pisa-)tower of states associated with the vac- 
uum of magnons |0) (these states are called QDJS for 
Quasi Degenerate Joint States in refs. ||,^), and the 
tower of states associated with the softest magnon |1). 
The semi-classical Neel ground-state with an order pa- 
rameter pointing in the {0, </>) direction is a wave packet 



of eigen-levels with different total spin belonging to the 
QDJS |0). Mixing different S values up to VTV allows the 
localization of the order parameter. If the levels of the |0) 
family collapse to the ground-sate faster than the softest 
magnon, the classical picture survives to quantum fluctu- 
ations and the system exhibits LRO: as for the solid, that 
might be the case at T = for 2d magnets. For T ^ 0, 
thermal fluctuations destroy the order (Mermin- Wagner 
theorem). Notice that the entropy of the semi-classical 
ground-state scales as Ln{N'^), and the entropy per spin 
of the ordered ground-state is zero in the thermodynamic 
limit. 

The levels associated to the vacuum of magnons have 
special properties with respect to the lattice symmetry 
group: these properties reveal the space symmetry break- 
ing of the Neel state. Let us take as an example the 
coUinear order on the honeycomb lattice (Fig. 3): it 
does not break the translational symmetry of the Bra- 
vais lattice (thus only the k = momentum appears in 
the QDJS). But it breaks both the inversion symmetry 
with respect to the center of an hexagon and the reflexion 
with respect to an axis going through nearest neighbor 
hexagon centers, and it is invariant in the product of 
these two symmetries. Thus both the even and odd ir- 
reducible representations of these two symmetry groups 
appear in the QDJS (symbols Rt; — ±1, a = ±1), but 
they are always associated so that 



R-^.o- 



-1 



(4) 



which insures the invariance in the product of the two 
symmetries, and 

a.(-l)^ = +l (5) 

which insures the invariance of the semi-classical ground- 
state in the product of a reflexion times a spin-flip. This 
is the simplest situation for an antiferromagnet. On the 
square lattice (which is equally bipartite) Neel order is 
coUinear and Eq. 3 remains valid, but as the translation 
symmetry is broken, both the k — (0, 0) and k — (tt, tt) 
vectors appear among the QDJS. 

A more complex order parameter implies a richer struc- 
ture of the QDJS: in the case of the Heisenberg model 
on the triangular lattice the order parameter is no longer 
a vector but a trihedron (or equivalently a rotation ma- 
trix). Eq.3 is then more complex, involving both the 
in-plane and out of plane susceptibilities, and there is 
an internal generation of an extra quantum number (the 
projection of the total spin on the helicity of the magnet). 
The free precession hamiltonian is that of a symmetric 
top m|,|ll| and 0{N^) states collapse to the ground-state 
in the thermodynamic limit. (Quantum mechanically the 
QDJS associated to p-sublattice Neel order are adiabat- 
ically connected to the states obtained by addition of p 
spins of size 



^This energy has to be measured along the arrow in Fig 1: 
this takes care of the momentum of the phonon. 



N_ 
2p- 



A thorough analysis of the symme- 
try breaking associated to Neel order on the triangular 
lattice is done in fBl . 



Unexpectedly the mechanism of "order by disorder" 
can be detected by looking at the finite size scaling of 
relatively small samples [|2| . We investigated this mech- 
anism in the Ji — J2 model on the triangular lattice. 
For 1/8 < J2/J1 < 1 the classical ground-state ex- 
hibits a large degeneracy associated to different kinds 
of 4-sublattice Neel order. For small sizes {N = 16) 
the quantum fluctuations are weak and the vacuum of 
magnons involve all the levels associated with addition 
of four momenta. In the N ~ 28 sample quantum fluc- 
tuations are already strong enough to lift the degeneracy 
of this family: the levels characteristic of the 2-sublattice 
Neel order remain the lowest ones in the spectrum, the 
other ones drift towards the continuum of excitations. It 
is a nice illustration of the fact that long wave length 
quantum fluctuations tend to restore the symmetries of 
the original hamiltonian. The quantum ground-state is 
more symmetric than the classical one. 

Order parameters, spin susceptibilities [^, stiffnesses 
|lO| can be extracted from the numerical data: they are 
interesting pieces of data to be compared to the spin 
waves calculations at finite sizes, but as it has been dis- 
cussed in the beginning of this section, the leading term of 
their finite size scaling is 0(A^^2 ) and therefore it is more 
significant to look to the tower of QDJS which collapse as 
0{N~^). It is also simpler to look at the symmetries of 
the low lying levels and deduce the symmetry of the order 
parameter than to determine it by direct computations. 
As a consequence when looking at the disappearance of 
Neel LRO due to extra frustration or quantum fluctua- 
tions it is preferable to pinpoint the range of parameters 
where the tower of states is destroyed than search for the 
vanishing of the order parameter ||l^;n . 



III. DIMER LONG RANGE ORDER AND 
VALENCE BOND CRYSTALS 

The second class of long-range ordered systems does 
not imply SU{2) symmetry breaking but only space sym- 
metry breaking. We will call these systems Valence Bond 
Crystals (VBC) g: their ground-states are characterized 
by regular arrays of singlet states on simple bonds or 
more complicated subsets (4-site plaquettes...). In the 
thermodynamic limit the ground-state of a VBC is de- 
generate, with a finite degeneracy directly related to the 



■^The name Valence Bond Solid is already "registered" by Af- 
fleck, Kennedy Lieb and Tasaki ||la,t6| for very specific sys- 
tems where the individual spins obey the relation 2S = z 
(with z the coordination number of the lattice). In such cir- 
cumstances there exists a whole class of Hamiltonians which 
have a unique highly symmetric ground-state where each site 
shares a singlet with any neighbor connected to it. 



symmetry breaking. The VBC on the Ji — J2 chain (for 
^ > 0.2411...) is doubly degenerate (it breaks the one- 
step translation symmetry of the lattice). The simplest 
VBCs that could be encountered on a triangular lattice 
have a 12-fold degeneracy (a factor 4 comes from the 
translation symmetry breaking and a factor 3 from the 
breaking of C3). It is a common belief that a VBC has a 
gap to all excitations (except indeed at critical points). 

In 2d, the more convincing evidence of such a phase, 
with a spontaneous discrete symmetry breaking has been 
seen on the Ji — J2 model on the square lattice (see [p|jl7|] 
and refs. therein). In that case the VBC breaks trans- 
lational symmetry but probably not C4: it is probably a 
long range order of S" = plaquettes. 

We have a new example of a VBC in the spectrum of 
the Ji — J2 model on the honeycomb lattice. When the 
second neighbor coupling is strong enough the frustration 
destroys the collinear Neel long range order and gives 
birth to a pattern of valence bonds (with a Z3 symmetry 
breaking associated to the three spatial directions for the 
dimerized bonds). Our predictions is based on the scahng 
of the low lying singlets of the spectra, it is in agreement 
with the dimer-dimer correlations functions measured in 
the first singlet levels supposed to be degenerate in the 
thermodynamic limit p]. 

The nature of the first excitations of these models (co- 
herent or incoherent magnons) remains to be elucidated. 



IV. FULLY GAPPED SPIN-LIQUIDS 

Besides these situations with symmetry breaking and 
LRO in a local observable, we have found two other 
behaviors, probably generic, with no obvious symmetry 
breaking, no long range order in a local order parameter. 
We call these situations "Spin-Liquids" and will succes- 
sively describe both of them. 

Let us first describe very quickly the simplest one 
which has the following properties (an extensive descrip- 
tion of the properties of this phase is given in |lq,n9|): 

• A gap to any excitations 

• All spin-spin correlations are short ranged 

• A four-fold degeneracy of the ground-state on a 2- 
torus 

We have observed this phase for the first time in the 
multiple-spin exchange problem. This effective hamil- 
tonian has been introduced by Eyring and Thouless to 
describe the magnetism of quasi localized fermions (elec- 
trons in the Wigner Crystal on one hand, ^He atoms 
in solid phase on the other). The Hamiltonian of the 
multiple-spin exchange model is given by 

H = ^(-1)" J„(P„ + P-I) , J„ > , n > 2 (6) 



where J„ are the n-spin exchange tunnehng probabilities 
(exchange coefficients), Pn and P~^ are the n-spin per- 
mutation operators and their inverse. The alternate sign 
in the summation over n in Eq. ra comes from the permu- 
tation of fermions. In general, the exchange coefficients 
decrease with increasing n. The two-spin exchange term 
gives the Heisenberg Hamiltonian up to a constant: 



Po — 2S,; • So 



(7) 



where S^ and Sj are spins localized at site i and j, re- 
spectively. The three-spin exchange operator is equal to 
a sum of two-spin exchange operators |19| . Thus, the 
two and three-spin exchange terms are described (up to 
a constant) by an Heisenberg hamiltonian with an ef- 
fective coupling constant Jf*^ = J2 — 2J3. Cyclic four- 
spin exchange gives birth to new physics: it is quartic in 
spin operators and has built-in frustration. As a result 
the ground-state of this 4-spin exchange hamiltonian on 
a single plaquette is doubly degenerate with total spin 
5 = and 1! At the classical level this hamiltonian 
displays various antiferromagnetic phases pO|-p^ , which 
are wiped out by quantum fluctuations, leading in a large 
range of parameters to the Spin-Liquid with the above- 
mentioned properties |18| , |l9|| . 

The four-fold degeneracy of the ground-state on a 2- 
torus is a topological property which was predicted more 
than ten years ago pi^-p5[|. Contrarily to Oshikawa ||2q ] 
we do not think that it is the signature of long range order 
in a local order parameter (see discussion in Misguich 

Very few things are clearly stated on the excitations of 
this Spin-Liquid: 

• Plateaus are expected in the magnetization curves 
both for zero and half the saturation magnetization. 

• From the finite size spectra it might be suspected 
that the triplet excitations form a "continuum" 
above the spin-gap and that an S* = bound-state 
exists just below the spin-gap. Up to now, a clear 
signature of deconfined spinous has not been found 
in this Spin-Liquid. Futher investigations appear 
necessary. 

As a further hint that might be of importance: short 
range correlations are slightly ferromagnetic. In fact 
this Spin-Liquid phase appears in the vicinity of a fer- 
romagnetic instability, as can be seen both in the case of 
the multiple-spin exchange Hamiltonian on the triangular 
lattice, but also for the Ji — J2 model on the honeycomb 
lattice ||]. 

Experimentally, this phase might be observed in the 
low density 2d layers of solid '^He where there is evidence 
that the coupling constants might be in the good range 
(see |19| , p7| -p9[ and references therein) and possibly in the 
Wigner crystal pO|. 



V. KAGOME-LIKE SPIN-LIQUIDS 

The last situation that we will now discuss is by some 
details more puzzling. It has been obtained by deplet- 
ing the triangular lattice, removing one spin out of four: 
the lattice which is obtained is called the kagomc lattice. 
Some possible physical realizations are SrCrGaO or var- 
ious jarosites |3^-|3^]. The model has been extensively 
studied in the classical limit [p4|-p8[ and it displays an 
extensive local degeneracy (i.e. an extensive number of 
modes of zero energy). Some authors have pleaded in 
favor of a selection of a particular ground-state by the 
"order out disorder" mechanism ^^1 but the finite size 
scaling of spin-1/2 spectra contradicts this assumption 
|13| . The system is certainly a Spin-Liquid in the follow- 
ing sense: 

• all the correlations (spin-spin, dimer-dimer, chiral- 
chiral, ...) are very small and decay with distance 
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• the system has a gap for magnetic excitations. 
Even if it is not sure that the largest size (N=36) 
is much larger than all correlation lengths Q , the 
finite size scaling studies point to a small spin-gap: 
probably of the order of one tenth of the coupling 
constant and quite likely larger than one twentieth 
of the coupling constant |l^,Q . 

• the spin-gap is filled by a continuum of singlets and 
the number of these singlets increases exponentially 
with the system size A^ PjiH ^^ shown in Fig. 4. 

The understanding of the low lying continuum of sin- 
glets is still incomplete but important progresses have 
been done recently |{43[-[45|] . Mambrini and Mila Q] have 
shown that the low lying singlet levels could be related to 
a family of short range dimer coverings of the lattice. It 
must be emphasized that short range dimer coverings on 
this lattice is not necessarily synonymous with exponen- 
tial decreasing correlations and gapful excitations: the 
recent results of Mambrini and Mila prove that this is 
indeed not the case and the comparison between these 
results and those of Zeng and Elser Q point to the fact 
that the non-orthogonality of the dimer coverings plays 
a major role in this respect. This is an important step 
forward and we would like now a simple picture of the 
first excitations. 

Ten years ago, Rokhsar and Kivelson have stud- 
ied a model of hard-core dimcrs on the square lat- 
tive p3| , p6| -^. The phase diagram of their model dis- 
plays two VBC phases separated by a critical point. At 



The spin-gap of the N = 21, 27, 36 samples seems to be in 
a cross-over regime of scaling, but the energy per spin still 
evolves as N~^^^\ 



the transition point between the two phases the spectrum 
of excitations becomes gapless. The low lying modes are 
the Goldstone modes associated to a quasi LRO in the 
orientation of the dimers. Let us assume that this hy- 
pothesis is qualitatively valid on the kagome lattice and 
derive its consequences for the number of states below 
the spin-gap. Assuming that these Goldstone modes are 
non-interacting, it is easy to extract from their disper- 
sion law {e{k) — fc"), the internal energy and entropy 
versus temperature and thus the number of states below 
the first triplet excitation. The logarithm of this number 
of states scales as N'^^+^ . In Fig. 4 we compare the nu- 
merical results with this predicted law (assuming n = 1 
or 2). Fixing the unknown constants from the measured 
value for A^ = 12, we see that the prescription is an order 
of magnitude off at TV = 36: the n — 1 mode description 
gives ^ 20 singlet levels in the spin-gap, whereas their 
real number is 210 [|! With the present pieces of infor- 
mation we must conclude that the description of the low 
lying singlets as Goldstone modes associated to a quasi 
LRO in the dimer covering is not supported by the nu- 
merics Q A simple picture of these first excitations is still 
lacking! 

On the other hand, the exponential increase in the 
number of eigenlevels in a finite band of energy implies 
a finite entropy per spin at T=0. The same phenomena 
appearing in each spin sector might be at the origin of the 
anomalous glassy behavior "without chemical disorder" 
seen in various compounds (see p9| and refs therein). 

Less speculative, the presence of a finite entropy in the 
singlet channel at low temperature is consistent with the 
observation of Ramirez and coworkers |E3| that the low 
temperature specific heat of SrCrGaO is essentially in- 
sensitive to large magnetic fields (much larger than the 
temperature) ||5l| . This picture is also compatible with 
the muons experiments of Uemura et al %% and the elas- 
tic spin diffusion measurement of Lee et al |53|1 . 

To conclude wc must add that these systems seem able 
to support unconfined spinous excitations with non zero 
Chern numbers |{42| . 

Up to very recently, we were thinking that this second 
kind of Spin-Liquid was associated to special lattices as 
the kagome or pyrochlore lattices, which at the classical 
level display an infinite local degeneracy. Some recent re- 



* This number is slightly larger than the one quoted in M] 
(some degeneracies were incorrectly numbered). 

^A second check has been done on the density of states ver- 
sus E for A^ = 36. Due to the intrinsic quantum noise of the 
numerical density of states, it is impossible with this last cri- 
terium to really discriminate between the different hypothe- 
ses, but the density of states implied by a mode description 

2 

exp{aE "+2 ) seems less probable than a pure exponential or 
power law ( the x^ criterium is approximately twice as large) . 



suits bring new pieces of information: we have observed 
kagome-like spectra on the triangular lattice when the 
3-sublattice Neel order is destabilized by a frustrating 
4-spin exchange p4| and we have no evidence of an infi- 
nite local degeneracy at this point. We now conjecture 
that the interesting feature common to these two cases 
is the appearance of the " kagome-like" phase after desta- 
bilization of a non-coUinear long-range order (either due 
to quantum fluctuations -when going from the triangu- 
lar to the kagome lattice by weakening of bonds u^-, or 
to competing frustrating interactions in the multiple-spin 
exchange problem |H). Contrarily, destabilization of a 
coUinear order (like in section II) leads to dimer LRO. 

Is this " kagome-like behavior" representative of a true 
new phase or of a quantum critical point? This remains 
an open question. In the two studied cases we found these 
characteristics in a rather wide range of parameters: It is 
nevertheless impossible to exclude a finite size effect. If 
criticality is to be viewed as a situation where there is no 
intrinsic length scale nor energy scale and a slowing down 
of all dynamics, it remains unclear if the proliferation of 
low lying singlets in this model could not be a signature 
of criticality (and the above assumption on the scaling of 
the density of states much too naive!). 



VI. MISCELLANEOUS REMARKS AND 
CONCLUSION 

As a last remark we would like to discuss some ques- 
tions around the words and concepts of order and disor- 
der. The established expression of " order by disorder" (or 
" order out of disorder" ) expresses the fact that quantum 
fluctuations (as well as thermal fluctuations for which it 
was primarily introduced by Villain) stabilize the simpler 
and the more symmetric ordered state amongst the al- 
lowed r = classical ground-states. But more generally 
the effect of long wave-length quantum fluctuations on 
a (semi-)classical ground-state might lead to still simpler 
and more symmetric situations which do not break SU{2) 
(section III) and even any space symmetries (sections IV 
and V). In such an issue the word "quantum disorder" 
has been used (by opposition to Neel long range order) 
but indeed it should be used with great caution: the 
ground-state of the Spin-Liquid of Section IV is highly 
symmetric and very far from a disordered system. Only 
the "kagome-like" state might be qualified as "quantum 
disordered" (section V). To measure this delicate issue 
one may look at the T — entropy per spin of the ground- 
state. In the thermodynamic limit, this quantity scales 
as O {pLn{N)/N) for a Neel ordered state with p sublat- 
tices, as 0{1/N) in the Valence Bond Crystals (Sec. Ill) 
and in the fully gapped Spin-Liquid (Sec. IV), and it is 
0(1) in the "kagome-like" Spin-Liquids. 

In conclusion exact diagonalizations have allowed a 
characterization of 4 different and (probably ) generic 



ground-states of quantum antiferromagnets in 2d. The 
two kinds of Spin-Liquids conjectured by P.W. Ander- 
son in 1973 Q have been found (the second with the 
continuum of singlets seems nevertheless rather different 
from the initial conjecture). Possible physical realiza- 
tions of these Spin-Liquids exist. There remain many 
questions on the first excitations of these Spin-Liquids. 
Some progress on this subject can probably be done in 
computing dynamical structure functions with exact di- 
agonalizations schemes, but it would also be useful to 
devise new analytical approximate tools to look at these 
questions: exact results on the ground-state and low ly- 
ing levels of the spectra may hopefully help in achieving 
this goal! 
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FIG. 1. Typical spectrum of a finite size solid. The first 
tower of eigenlevels |0) is associated with c.o.m. motion. The 
second tower |1) describes the softest phonon. 
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FIG. 3. Antiferromagnetic Heisenberg model on the hon- 
eycomb lattice: eigenenergies vs the eigenvalues of S^. The 
dashed-line is a guide to the eyes for the QDJS. The dotted 
line joins the states associated to the first magnon. There is 
one QDJS for each S (as expected for a coUinear antiferromag- 
net): they are k = states, invariant under a 27r/3 rotation 
around an hexagon center, even (odd) under inversion, odd 
(even) under a reflexion with respect to an axis going through 
nearest neighbor hexagon centers for S even (odd). 
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FIG. 2. Typical spectrum of a finite size coUinear anti- 
ferromagnet with a Neel order. The tower of eigenlevels |0) 
are the QDJS associated with the free dynamics of the order 
parameter. The second tower |1) is associated with the lowest 
magnon. 
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FIG. 4. Logarithm of the number of states in the mag- 
netic gap on the kagome lattice (squares fitted by a continous 
line) vs the size A*' of the system. The observed behavior 
differs from the scaling law exp{bN"'^"^^^) deduced from a 
single mode description of the continuum with a dispersion 
law e{k) — k" for n = 1 (short-dashed line) and n — 2 
(long-dashed line). 



